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Analyzing correlation functions of charmonia at finite temperature (T) on 32 3 x(32— 96) anisotropic lattices 
by the maximum entropy method (MEM), we find that J/ip and r/c survive as distinct resonances in the plasma 
even up to T ~ 1.6T C and that they eventually dissociate between 1.6T C and 1.9T C (T c is the critical temperature 
of deconfinement). This suggests that the deconfined plasma is non-perturbative enough to hold heavy-quark 
bound states. The importance of having sufficient number of temporal data points in the MEM analysis is also 
emphasized. 



1. INTRODUCTION 

Whether hadrons survive even in the decon- 
fined quark-gluon plasma is one of the key ques- 
tions in quantum chromodynamics (QCD). This 
problem was first examined in pp and [2] in differ- 
ent contexts. The fate of the heavy mesons such 
as J/tjj in the deconfined plasma was also inves- 
tigated in a phenomenological potential picture 
taking into account the Debye screening [J . In 
general, there is no a priori reason to believe that 
dissociation of the bound states should take place 
exactly at the phase transition point 4 . 

From the theoretical point of view, the spectral 
function (SPF) at finite temperature T, which has 
all the information of in-medium hadron proper- 
ties, is a key quantity to be studied. Recently, 
the present authors have shown J5J that the first- 
principle lattice QCD simulation of SPFs is pos- 
sible by utilizing the maximum entropy method 
(MEM). We have also formulated the basic con- 
cepts and applications of MEM on the lattice at 
T = and T ^ in 6 . To draw the conclu- 
sion with a firm ground, in the following we put 
special emphases on [I] the MEM error analysis 
of the resultant SPFs and [II] the sensitivity of 
the SPFs to A^ata (the number of temporal data 
points adopted in MEM). These tests are crucial 
to prevent fake generation and/or smearing of the 
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peaks and must be always carried out as empha- 
sized in [51715] , 



2. LATTICE SIMULATION 

The applications of MEM to T ^ system have 
been known to be a big challenge HJ. The diffi- 
culty originates from the fact that the temporal 
lattice size L T is restricted as L T = l/T = N T a T , 
where a T (N T ) is the temporal lattice spacing 
(number of temporal lattice sites). Because of 
this, it becomes more difficult to keep enough 
A'data to obtain reliable SPFs as T increases. 
Thus, simulations up to a few times T c with A^ata 
as large as 30 |7J inevitably require an anisotropic 
lattice. 

On the basis of the above observation, we 
have carried out quenched simulations with [3 — 
7.0 on 32 3 x N T anisotropic lattice with N T = 
32, 40, 46, 54, and 96. The renormalized 
anisotropy is £ = a a /a T = 4.0 with a a being the 
spatial lattice spacing. We take the naive plaque- 
tte gauge action and the standard Wilson quark 
action. Corresponding bare anisotropy £o = 3.5 
is determined from the data given in |9I10| . The 
fermion anisotropy = n T /n a = 3.476 with 
n a = 0.08285 (k t = 0.2880) being the the spa- 
tial (temporal) hopping parameter is determined 
by comparing the temporal and spatial effective 
masses of pseudoscalar and vector mesons on a 
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32 2 x 48 x 128 lattice. a T = a a /A = 9.75 x 1(T 3 
fm is determined from the p meson mass in the 
chiral limit. The masses determined on T = 
(32 2 x 48 x 128) lattice are m J/4 , = 3.10 GeV and 
m rjc = 3.03 GeV. If the temporal distance be- 
tween the source and the sink is closer than £a T , 
lattice artifact due to anisotropy would appear in 
the SPFs for co > ir/^a T . To avoid this, we ex- 
clude the six points near the edge (n = 1,2,3 
and N T — 3, N T — 2, N T — 1) and adopt the points 
Ti = 4, 5, • • • and N T — 4, N T — 5, • • • until we reach 
the total number of points Adata(< N T — 7) . 

3. RESULTS 

Shown in Fig^are p{w)s 2 for J/ip at T/T c = 
0.78,1.38, and 1.62 (FigDJa)) and those at 
T/T c = 1.87 and 2.33 (FigQ3») (Correspond- 
ing A^data used in these figures are 7V data = 
89, 40, 34, 33, and 25 from low T to high T). 
If the deconfined plasma were composed of al- 
most free quarks and gluons, SPFs would show 
a smooth structure with no pronounced peaks 
above the qq threshold. To the contrary, we find a 
sharp peak near the zero temperature mass even 
up to T ~ 1.6T C as shown in Fig^a), while 
the peak disappears at T ~ 1.9T C as shown in 
FigGJb). The width of the first peak in Fig^a) 
partly reflects the unphysical broadening due to 
the statistics of the lattice data and partly reflects 
possible physical broadening at finite T. At the 
moment, the former width of a few hundred MeV 
seems to dominate and we are not able to draw 
definite conclusions on the thermal mass shift and 
broadening. 

Let us now evaluate the reliability of the ex- 
istence (absence) of the sharp peak at T/T c = 
1.62(1.87) by the two tests [I] and [II] mentioned 
before. First test is the error analysis of the 
peak. Shown in Fig [5] are the SPFs for J/ip at 
T/T c = 1.62 and 1.87 with MEM error bars (The 
frequency interval over which the SPF is averaged 
is characterized by the horizontal position and ex- 
tension of the bars, while the mean value and the 
la uncertainty of the integrated strength within 
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-^or further details of the lattice simulation, see 151 . 
2 Following |5J, we define dimensionless SPFs: A(uj) 
ui 2 p(w) for rjc and A(a>) = 3io 2 p(io) for J/t/j. 



Figure 1. Spectral functions for J/ip (a) for 
T/T c = 0.78,1.38, and 1.62 (b) for T/T c = 1.87 
and 2.33. 



the interval are characterized by the heights of 
the bars). The sharp peak at T — 1.62T C is sta- 
tistically significant, and the absence of the peak 
at the same position at T = 1.87T C is also statis- 
tically significant. The features found for J/ip in 
Figs. n an d0are also observed for r\ c . 

The second test is the A^ata dependence of 
the SPFs. Shown in Fig|3Ja) is the comparison 
of SPF obtained with Adata = 34 and that of 
A^data = 39 for the same temperature T = 1.62T C 
(N T = 46). The two curves are almost identi- 
cal with each other. Fig|3Jb) shows SPFs ob- 
tained with A^data = 26 and 33 for higher tem- 
perature T = 1.87T C (AV = 40). Again the two 
curves are almost identical. Therefore, the quali- 
tative change of the SPF between T — 1.62T C and 
1.87T C is the real thermal effect and not caused 
by the artifact of the insufficient number of data 
points. 

To conclude, we have studied the spectral func- 
tions of J/ip and r\ c in the deconfined plasma us- 
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Figure 2. Spectral functions for J/ip with MEM 
errors (a) for T/T c = 1.62 (b) T/T c = 1.87. 



ing lattice Monte Carlo data and the maximum 
entropy method. The number of temporal sites 
N T is taken as large as 46 and 40 for T/T c = 1.62 
and 1.87, respectively. Careful analyses of the 
MEM errors and the A^ata dependence of the re- 
sults are carried out. It is found that there are 
distinct resonances up to T ~ 1.6T C and they dis- 
appear between 1.6T C and 1.9T C . 
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